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Abstract  The  rendezvous  and  formation  problem  is  a  significant  part  for  the  unmanned  aerial 
vehicle  (UAV)  autonomous  aerial  refueling  (AAR)  technique.  It  can  be  divided  into  two  major 
phases:  the  long-range  guidance  phase  and  the  formation  phase.  In  this  paper,  an  iterative 
computation  guidance  law  (ICGL)  is  proposed  to  compute  a  series  of  state  variables  to  get  the  solu¬ 
tion  of  a  control  variable  for  a  UAV  conducting  rendezvous  with  a  tanker  in  AAR.  The  proposed 
method  can  make  the  control  variable  converge  to  zero  when  the  tanker  and  the  UAV  receiver  come 
to  a  formation  flight  eventually.  For  the  long-range  guidance  phase,  the  ICGL  divides  it  into  two 
sub-phases:  the  correction  sub-phase  and  the  guidance  sub-phase.  The  two  sub-phases  share  the 
same  iterative  process.  As  for  the  formation  phase,  a  velocity  coordinate  system  is  created  by  which 
control  accelerations  are  designed  to  make  the  speed  of  the  UAV  consistent  with  that  of  the  tanker. 
The  simulation  results  demonstrate  that  the  proposed  ICGL  is  effective  and  robust  against  wind 
disturbance. 
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1.  Introduction 


With  the  rapid  development  of  UAV  technologies,  UAVs  are 
being  used  to  carry  out  various  missions  like  high  altitude 
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surveillance  and  reconnaissance,  long  distance  military  strikes, 
etc.,  some  of  which  require  UAVs  to  continuously  stay  in  the 
air  for  a  quite  long  time.  To  this  end,  autonomous  aerial  refu¬ 
eling  (AAR)  becomes  a  key  issue  to  be  addressed  for  these 
applications  of  UAV.  There  are  four  main  tasks  to  be  accom¬ 
plished  sequentially  during  AAR  of  UAV:  UAV  rendezvous, 
formation  maintaining,  pipeline  docking,  and  refueling.  This 
paper  mainly  focuses  on  presenting  a  guidance  law  for  the  ren¬ 
dezvous  and  formation  maintaining  processes.  For  decades,  a 
lot  of  achievements  have  been  made  to  develop  guidance  tech¬ 
nologies  for  the  UAV  rendezvous  problem.  One  of  the  most 
mature  technologies  is  the  classical  proportional  navigation 
guidance  (PNG). 4  After  that,  many  improved  proportional 
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navigation  guidance  (IPNG)  schemes  have  been  evolved  based 
on  PNG.  To  further  enhance  the  guidance  performance, 
modern  control  theories,  such  as  sliding  mode  control,8  dif¬ 
ferential  geometric  method,11  13  neural  network,14  lo  and  so 
on,  are  also  employed  for  the  design  of  UAV  rendezvous  guid¬ 
ance  laws.  For  formation  maintaining,  there  are  several  exist¬ 
ing  methods.  In  Ref.  17,  a  full-state  linearization  via  a  dynamic 
feedback  controller  is  designed  for  controlling  two  robots  in  a 
leader-follower  configuration.  In  Ref.  18,  a  synchronized  posi¬ 
tion  tracking  controller  is  incorporated  in  formation  flight  con¬ 
trol  for  multiple  flying  wings.  Ref.  19  proposed  a  new 
approach  of  hybrid  supervisory  control  for  the  leader-follower 
formation  problem.  The  hybrid  supervisory  control  approach 
provides  a  tractable  framework  for  hybrid  synthesis  of  forma¬ 
tion  control.  Within  this  framework,  a  new  method  of  abstrac¬ 
tion  based  on  polar  partitioning  of  the  state  space  is 
introduced.  Ref.  20  presented  an  iterative  guidance  method 
for  launch  vehicles.  In  this  method,  the  guidance  for  a  launch 
vehicle  is  formulated  as  an  optimal  control  problem,  in  which 
the  transient  state  of  the  vehicle  is  taken  as  the  initial  value  and 
the  target  point  as  the  terminal  constraint.  The  objective  func¬ 
tion  is  to  minimize  the  flight  time  of  the  vehicle  moving  from 
the  current  position  to  the  target  point.  During  the  whole  flight 
process  of  the  vehicle,  for  each  time  interval,  the  control  solu¬ 
tion  and  the  corresponding  flying  trajectory  are  obtained  by 
solving  the  established  guidance  equations.  Through  a 
repeated  iterative  computation,  the  launch  vehicle  is  eventually 
guided  to  the  target  point  and  satisfies  its  predefined  state.  This 
method  has  been  studied  for  the  guidance  and  control  design 
of  launch  vehicles  and  ballistic  missiles.21  24  However,  these 
approaches  involve  heavy  computation  load  and  are  difficult 
for  practical  engineering  applications. 

In  this  paper,  a  novel  iterative  computation  guidance  law 
(ICGL)  is  proposed  for  a  UAV  to  perform  rendezvous  and  for¬ 
mation  maintaining  with  a  tanker  in  the  AAR  process.  In  the 
ICGL,  the  rendezvous  process  of  a  UAV  with  a  tanker  in 
AAR  is  divided  into  two  major  phases:  the  long-range  guidance 
phase  and  the  formation  phase.  In  each  phase,  the  ICGL  com¬ 
putes  the  relative  state  parameters  between  the  UAV  and  the 
tanker  in  real-time  to  obtain  the  specific  control  variable.  Fur¬ 
thermore,  the  ICGL  divides  the  long-range  guidance  phase  into 
two  sub-phases  which  share  the  same  iterative  process:  the  cor¬ 
rection  sub-phase  and  the  guidance  sub-phase.  The  ICGL  pro¬ 
posed  in  this  paper  is  totally  different  from  the  above 
mentioned  iterative  guidance  method  based  on  the  optimal  con¬ 
trol  theory.  In  the  ICGL,  a  UAV  approaches  to  a  tanker  along 
a  smooth  arc  trajectory  designed  by  a  geometric  method.  Thus 
the  iteration  is  based  on  the  relative  position  between  the  UAV 
and  the  tanker  in  each  time  interval.  The  designed  arc  trajectory 
strategy  has  an  advantage  which  makes  the  control  vector  per¬ 
pendicular  to  the  velocity  of  the  UAV.  Therefore,  it  is  easy  to  be 
realized  in  practical  engineering  perspective. 

The  remainder  of  the  paper  is  organized  as  follows.  In  Sec¬ 
tion  2,  the  problem  of  rendezvous  for  AAR  is  formulated 
mathematically.  In  Section  3,  the  ICGL  algorithm  is  developed 
for  a  UAV  to  perform  rendezvous  and  formation  maintaining 
with  a  tanker  in  AAR.  In  Section  4,  simulations  are  performed 
to  verify  the  effectiveness  of  the  proposed  ICGL  by  comparing 
it  with  the  nonlinear  guidance  (NG)  method.  Then  simulation 
of  wind  disturbance  injection  is  conducted  to  demonstrate  the 
robustness  of  the  ICGL.  The  conclusion  remarks  are  given  in 
Section  5. 


2.  Problem  description 

An  illustrative  diagram  for  the  rendezvous  process  of  a  UAV 
with  a  tanker  for  aerial  refueling  is  given  in  Fig.  1.  The 
research  objective  of  this  paper  is  to  guide  the  UAV  near  the 
tanker  and  then  the  UAV  keeps  a  formation  flight  with 
the  tanker. 

In  the  UAV  guidance  process,  there  exist  two  situations  for 
the  relationship  of  the  initial  velocity  between  the  UAV  and 
the  formation  location:  the  unparallel  one  and  the  parallel 
one.  The  system  can  be  described  in  the  Earth-fixed  inertial 
frame  Ocxyz  representing  a  general  unparallel  situation  shown 
in  Fig.  2  with  its  origin  at  the  initial  location  of  the  UAV 
denoted  as  Oc.  The  orientation  of  the  x-axis  is  consistent  with 
the  initial  velocity  vector  of  the  UAV.  The  z-axis  is  downward 
and  perpendicular  to  the  x-axis  in  the  vertical  plane.  The  y-axis 
is  determined  according  to  the  right-hand  rule.  Point  F  stands 
for  the  formation  location  and  point  P  is  the  projection  point 
of  F  on  the  Ocxy  plane.  H  is  an  intersection  point  at  which  PH 
is  perpendicular  to  the  x-axis.  The  transitional  state  of  the 
UAV  during  the  rendezvous  process  is  shown  in  Fig.  3,  in 
which  m  is  the  surface  which  passes  the  velocity  Vv  of  the 
UAV  and  is  parallel  with  the  x-axis  when  the  velocity  of  the 
UAV  is  unparallel  with  the  Ocxy  plane. 

The  position  vector  of  the  UAV  is  denoted  as 

P\j  =  Iau  5  Vu  5  ^u]T 

The  velocity  vector  of  the  UAV  is 

Vv  =  [vuxj  Vu y,  VUz ] 

The  motion  of  the  UAV  can  be  described  by  the  following 
kinematics  equations24 

(Xu  =  Vu  cos  G 

ju  =  vu  sin  a  cosy  (1) 

Zu  =  Vu  sin  g  sin  y 

where  Vu  is  the  magnitude  of  Vu  and 

VU  =  yj  tfjx  +  Vu y  +  v2Vz  (2) 

g  is  the  angle  between  the  velocity  of  the  UAV  and  the 
x-axis.  y  is  the  angle  between  the  surface  m  and  the  Ocxy  plane 
in  the  frame  Ocxyz  as  shown  in  Fig.  3. 

In  Fig.  2,  Fis  the  tracking  point  for  the  UAV  to  keep  a  for¬ 
mation  flight  with  the  tanker.  The  location  of  F  can  be  set 
according  to  the  current  position  of  the  tanker  with  a  deviation 
value  Ad  as  shown  in  Fig.  2.  The  position  vector  of  F  can  be 
expressed  as 


Fig.  1  An  illustrative  rendezvous  process  of  a  UAV  with  a 
tanker. 
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Fig.  2  Unparallel  situation. 


Fig.  3  State  of  the  UAV  during  the  guidance  process. 


P(=  P,~  \d=  [x{,y[,Z(]T 
The  velocity  vector  of  F  is 
Vi  =  Ft  =  [vfx,Vf,,,Vf.]T 

where  Pt  is  the  position  vector  of  the  tanker  and  Vt  is  the 
velocity  vector  for  the  tanker. 

The  motion  of  F  can  be  described  by  the  following 
kinematics 

{*f  =  Vfic 

yt  =  Vf v  (3) 

if  =  vfz 


3.  Design  of  guidance  law 

Since  the  parallel  situation  is  an  exceptional  case  of  the  unpar¬ 
allel  one,  the  unparallel  situation  is  mainly  discussed  in  this 
paper.  The  rendezvous  process  can  be  divided  into  two  phases: 
the  long-range  guidance  phase  and  the  formation  phase. 

3.L  Long-range  guidance  phase 

The  long-range  guidance  phase  is  defined  from  the  initial 
instant  to  the  instant  when  the  distance  between  the  UAV 
and  F  is  equal  to  or  smaller  than  a  given  value.  When  Vu 
and  Vf  are  not  parallel,  the  long-range  guidance  phase  should 
also  be  further  divided  into  two  sub-phases:  the  correction 
sub-phase  and  the  guidance  sub-phase. 


Fig.  5  Geometric  relationship  in  the  two-dimensional  plane. 


In  Fig.  5,  r  is  the  distance  between  the  UAV  and  P;  a  is  the 
angle  of  line-of-sight  from  the  UAV  to  P;  rj  is  the  advance 
angle  of  the  UAV;  Vp  is  the  velocity  of  projection  point  P 
on  the  Ocxy  plane;  i/j  is  the  angle  between  the  velocities  of 
the  UAV  and  P;  ex  is  the  unit  vector  of  the  v-axis  and  er  is 
the  unit  vector  along  the  direction  from  the  UAV  pointing 
to  P.  The  strategy  used  in  this  paper  to  guide  the  UAV  from 
the  initial  position  to  location  P  is  along  an  arc  orbital  approx¬ 
imation.  The  UAV  can  move  along  an  arc  trajectory  from  its 
current  position  to  the  target  location  P  by  applying  a  normal 
correction  overload  on  the  UAV. 

The  angle  a  of  the  UAV  satisfies  the  constraint 


a  = 


vu 

R 


where  the  radius  of  arc  trajectory  for  the  UAV  moving  to  P 
can  be  calculated  by 


2  sin  rj 

ex  •  ev 

where  n  =  a  —  a  ,  in  which  a  =  arccos  r—r. — r . 

kxlkrl 

According  to  the  analysis  above,  the  control  acceleration 
a\j  can  be  calculated  by 


(4) 


The  iterative  process  throughout  the  correction  sub-phase 
is  shown  below  in  detail. 


3.1.1.  Correction  sub-phase 

The  correction  sub-phase  is  to  adjust  the  direction  of  Vu  to  be 
parallel  with  the  velocity  of  the  projection  point  P  on  the  Ocxy 
plane  firstly.  The  parallel  situation  is  given  in  Fig.  4.  The  geo¬ 
metric  relationship  between  the  UAV  and  P  in  the  Ocxy  plane 
is  shown  in  Fig.  5. 


(1)  Initialization 

Initialize  the  following  parameters  at  the  initial  instant  t0 : 
The  position  vector  of  the  UAV 

PuM  =  [wj(M5ku(/o)>zu(/o)]T  =  [0,0, 0]T 
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The  velocity  vector  of  the  UAY 

VuM  =  [vUx(/o),  V\jy(to),  VUz(70)]  =  [Vu,0,0]T 

The  position  vector  of  P 

PpM  =  [xp(;0),jp(fo),Zp(fo)]T  =  [xp,>>p,0]T 
The  velocity  vector  of  P 
Vp(t0)  =  [Vpx(?o),  Vp y{t0),  Vpz(?o)]T  =  [vpx,  V.  v]T 

The  angle  o  of  the  UAV 

a(t0)  =  0 


The  relative  position  vector  of  P  to  the  UAY 
Los(/o)  =  Pp(*o)  —  Pu(t0) 

The  magnitude  of  the  distance  between  the  UAY  and  P 

r(*o)  =  |LosOo)l 

The  angle  of  line-of-sight 

/  \  Los(^o) 

a(*o)  =  arccos- — - — 

\ex\\r\os(t0)\ 

The  advance  angle  of  the  UAY 

rj(t0)  =  oc(t0)  -  a(to) 

The  magnitude  of  the  normal  correction  acceleration 

vu 2v^sin^0) 


av(to)  — 


R(k) 


r(to) 


(2)  Iterative  computation 

In  order  to  get  a  series  of  solutions  for  the  control  variable 
a\j  using  the  iterative  method,  the  differential  equations  in  Eqs. 
(1)  and  (3)  should  be  discretized  by  using  the  following 
formula 

...  dw  wit  +  A^)  —  wit) 

’'■«=W  A, 

By  this  way,  the  guidance  process  can  be  divided  into  multi¬ 
ple  discrete  sub-processes  as  shown  in  Fig.  6.  Each  sub-process 
is  regarded  as  a  straight  line  motion  with  time  interval  At.  As  is 
the  displacement  of  the  UAY  within  interval  At.  A is  the  cor¬ 
rection  angle  of  the  UAV  velocity  after  interval  At. 

From  the  discrete  sub-processes,  the  value  at  each  instant 
can  be  calculated  step  by  step.  Using  i  (i  =  1,  2,  3. . .)  represent 
the  ith  iterative  step,  the  iterative  computation  procedure  can 
be  illustrated  as  follows. 

Step  1.  Set  i  =  1  and  y  =  0. 

Step  2.  Calculate  the  position  vector  of  the  UAY  at  tt 
instant 

Pvfa)  = 


Fig.  6  Divided  guidance  process. 


where 


*u(A)  —  -^u(Ui)  +  v\jx(ti~i)At 


LufA')  —  LufA-i)  +  V\jy{ti-\)At 


\  zu (ti)  —  z\j(ti-\)  +  v\jz(ti-i)At 

Step  3.  Calculate  the  position  vector  of  P  at  tj  instant 

Pp(ti)  =  k(*).*('0,*pfo)]T 

where 


!Xp  (ti)  =  Xp  (f,_i)  +  Vpx(f;-i)A? 

yP(u)  =  +  vpv(;,_i)Af 

zP  {ti)  =  Zp0;-i)  +  vpz(ti-  \)\t 
Step  4.  Update  the  velocity  vector  of  P  at  tt  instant 

Vp  (U)  =  Vp(t,-i) 

Step  5.  Calculate  the  relative  position  vector  of  P  to  the 
UAV  at  ti  instant 

Los  (A')  =  Pp(ti)  —  P\]{tj) 


Step  6.  Calculate  the  magnitude  of  the  relative  position  vec¬ 
tor  of  P  to  the  UAY  at  tt  instant 

r(U)  =  |rlos(U)l 

Step  7.  Calculate  the  angle  of  line-of-sight  a  at  tt  instant 

/  \  Ct  Los  (t/) 

oc(ti)  =  arccos  - — n - --r 

\ex\\rios(ti)\ 

Step  8.  Update  the  angle  g  of  the  UAY  at  tt  instant 

c  (ti)  =  cr(u_i)  +  Afi(ti-i) 


Step  9.  Calculate  the  velocity  vector  of  the  UAV  at  tt  instant 

Vu(ti)  =  [Vux(ti)?  vUy(ti),  Vuz(ti)] 


where 


'  Vux(^)  =  Vu  COS  (7(A) 
<  Vuy(ti)  =  vusinff(t/) 

,  VUz(^)  =  0 


Step  10.  Calculate  the  advance  angle  of  the  UAV  at  tf 
instant 


n(fi)  =  «(U)  -  (7 (ti) 

Step  11.  Calculate  the  magnitude  of  the  UAV  normal  cor¬ 
rection  acceleration  at  tt  instant 

vu  _2v2,sin ti(ti) 
u(i)  R(t,)  r{tt) 


Step  12.  Calculate  the  correction  angle  of  the  UAY  velocity 
at  ti  instant 


mu)  = 


au^Af 

Vu 


Step  13.  Calculate  the  angle  between  the  velocity  vectors  of 
the  UAV  and  P  at  tt  instant 


[//( ti )  =  arccos 


Vu  (tt)  •  Fp  (ti) 

\vvmvvm 
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Step  14.  If  is  smaller  than  a  given  value  of  0.0005  rad 
in  this  paper,  terminate  the  iterative  computation  and  go  into 
the  guidance  sub-phase.  Otherwise,  let  i  =  i  +  1,  and  go  to 
step  2. 


3.1.2.  Guidance  sub-phase 

At  the  end  of  the  correction  sub-phase,  the  guidance  sub-phase 
is  activated.  At  the  same  instant  the  position  relationship 
between  the  UAV  and  F  as  well  as  P  in  the  frame  Ocxyz  is 
shown  in  Fig.  7. 

The  position  vector  of  F  in  Fig.  7  is 

Pf=  [xf,y[,Z(]T 

where  Xf  =  xp,  yf  =  yp,  zf  =  zd,  zd  is  the  required  position 
component. 

The  velocity  vector  of  F 

Vf  =  Vp 

In  order  to  use  the  iterative  computation  algorithm  similar 
to  that  used  in  the  correction  sub-phase,  a  coordinate  system 
Ogx'y'z'  is  created  with  its  origin  at  the  current  location  of 
the  UAV.  In  the  frame  Ogx'y'z',  the  orientation  of  the  V-axis 
is  consistent  with  the  current  velocity  direction  of  the  UAV. 
The  z'-axis  is  downward  and  perpendicular  to  the  V-axis  in 
the  vertical  plane.  The  /-axis  can  be  determined  according 
to  the  right-hand  rule  as  shown  in  Fig.  7.  The  coordinate  trans¬ 
formation  from  the  frame  Ocxyz  to  the  frame  Ogx'y'z'  can  be 
obtained  by 

P'  =  RgcP  +  D  (5) 


where  P  is  a  3  x  1  position  vector  in  the  Ocxyz,  P'  is  the  trans¬ 
formed  position  vector  in  the  Ogx'y'z'.  Rgc  is  the  rotation 
matrix  from  the  frame  Ocxyz  to  the  frame  Ogx'y'z'  which  can 
be  calculated  as  follows: 


R 


gc 


COS  (7 og 

—  sin  <Tog 
0 


sin  <70g 
COS  (J0g 

0 


0 

0 

1 


(6) 


where  <7og  is  the  value  of  o  when  the  frame  Ogx'y'z'  is  created 
and  D  is  a  3  x  1  displacement  vector: 

D  =  —RgcPog  (7) 

where  Pog  is  the  position  of  the  origin  Og  of  the  frame  Ogx'y'z' 
in  the  frame  Ocxyz.  Fig.  8  indicates  the  new  relationship 
between  the  UAV  and  F  in  the  frame  Ogx'y'z'.  For  the  guid¬ 
ance  sub-phase,  it  has. 

The  velocity  vector  of  the  UAV  in  the  frame  Ogx'y'z' 


Fig.  7  Position  relationship  at  the  end  of  the  correction  sub¬ 
phase. 


The  position  vector  of  F  in  the  frame  Ogx'y'z' 

P'f  =  RecP[  +  D=[X't,y'(,Z'c]T 

The  velocity  vector  of  F  in  the  frame  Ogx'y'z' 

The  iterative  computation  process  for  the  guidance  sub¬ 
phase  is  the  same  as  that  described  in  the  correction  sub-phase. 

3.2.  Formation  phase 

The  formation  phase  is  to  consider  how  to  keep  the  UAV  and 
the  tanker  flying  in  a  formation  for  the  preparation  of  the 
docking  of  the  refueling  probe  of  the  UAV  with  the  drogue 
of  the  tanker.  When  the  distance  between  the  UAV  and  the 
tanker  is  less  than  a  specified  value,  a  velocity  coordinate  sys¬ 
tem  Ovx''y"z"  is  created  as  shown  in  Fig.  9.  The  frame 
Owx"y"z"  is  used  to  measure  the  target  location  for  computing 
the  control  variable.  Fig.  9  shows  the  position  relation  between 
the  UAV  and  F  in  the  frame  Owx"y"z"  during  the  formation 
phase.  In  order  to  guide  the  UAV  to  the  formation  point  F 
and  keep  the  same  flight  direction,  a  control  acceleration  an 
is  designed  to  apply  on  the  UAV  with  direction  pointing  to 
Fi\  which  is  a  projection  of  formation  point  F  on  the  Owy"z" 
plane  in  the  velocity  coordinate  system  as  shown  in  Fig.  10. 
The  magnitude  of  an  is  determined  according  to  the  coordinate 
value  of  F  in  the  frame  Owx"y"z"  as 

an  =  A  sin  3  (8) 

where  A  is  a  constant  and  3  is  the  angle  between  the  UAV 
velocity  and  the  relative  position  vector  of  F  to  the  UAV. 

Besides,  the  formation  flight  phase  in  AAR  requires  the 
speed  of  the  UAV  to  be  consistent  with  that  of  the  tanker. 
For  this  reason,  the  control  acceleration  ah  (see  in  Fig.  10) 
along  with  the  direction  of  the  UAV  speed  should  be  utilized 
to  control  the  speed  of  the  UAV.  The  magnitude  of  ah  depends 


Fig.  8  New  position  relationship  in  Ogx'y'z'. 


Fig.  9  UAV  position  relative  to  the  formation  point. 
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Fig.  10  Control  direction  of  UAV. 


on  the  relative  distance  between  the  UAV  and  F  as  well  as  their 
current  speeds. 

The  flight  time  needed  for  the  UAV  to  cover  a  distance  of  s 
can  be  calculated  as 

2s 

t  = - 

Vu  +  Vf 

where  5  is  the  arc  trajectory  for  the  UAV  moving  to  F  com¬ 
puted  in  the  frame  Owx"y"z".  The  magnitude  of  the  control 
acceleration  ah  can  be  calculated  by 


In  summary,  the  rendezvous  process  for  the  UAV  with  the 
tanker  can  be  described  by  the  flow  chart  shown  in  Fig.  11. 

4.  Simulations  and  analysis 

In  order  to  verify  the  effectiveness  of  the  ICGL  for  the  rendez¬ 
vous  and  formation  maintaining  problem  in  AAR,  a  simula¬ 
tion  is  performed  here  using  the  proposed  ICGL  by 


comparing  with  the  NG  approach  presented  in  Ref.  25  to  dem¬ 
onstrate  the  effectiveness  of  the  ICGL. 

In  this  simulation,  assume  that  the  initial  conditions  of  the 
UAV  and  the  tanker  are  as  follows:Position  vector  of  the  UAV 

PuM  =  [*u(*o)3  Tu(*o)>  zv(k)]J  =  [0,0,0]Tm 
Velocity  vector  of  the  UAV 

VvM  =  [vux(^o)5  Vu-j,(^o),  Vuz(*o)]  =  [100,0,0]T  m/s 
Position  vector  of  the  tanker 

Pt(to)  =  [xt(t0),yt(e0),zt(e0)]T  =  [5000,  -3000, -1000]T  m 

Velocity  vector  of  the  tanker 

Vt(to)  =  [vtx(?o),vt,Vo),v,.(fo)]T  =  [50,  —40, 0]T  m/s 

During  the  formation  phase,  the  relative  position  of  the 
UAV  to  the  tanker  is  set  as 

A d=  [Adx,Ady,Adz]J  =  [35,  35,  50]T  m 

Simulations  for  the  rendezvous  and  formation  maintaining 
process  of  the  UAV  with  the  tanker  are  performed  by  using  the 
ICGL  and  the  NG  methods,  respectively.  Fig.  12  shows  the 
trajectories  of  the  UAV  approaching  to  the  tanker  in  the  ren¬ 
dezvous  and  formation  maintaining  process  using  the  ICGL 
and  the  NG. 

The  control  accelerations  for  the  correction  and  guidance 
sub-phases  are  shown  in  Fig.  13.  ac  and  ag  is  the  acceleration 
in  the  correction  sub-phase  and  guidance  sub-phase,  respec¬ 
tively.  Using  E  =  J 2  a2dt,26  we  can  compute  the  energy  con¬ 
sumptions  F^icgl  =  1108.5  for  the  ICGL  and  FNG  =  1270 
for  the  NG,  respectively.  Obviously,  the  UAV  consumes  less 
energy  to  accomplish  the  rendezvous  process  by  using  the 
ICGL  than  by  using  the  NG. 

The  altitude  change  histories  of  the  UAV  during  the  ren¬ 
dezvous  process  are  shown  in  Fig.  14  for  the  ICGL  and  the 


Fig.  11  Flow  chart  for  the  ICGL  approach. 
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Tanker  trajectory  - Using  ICGL - Using  NG 


Fig.  12  Trajectories  through  the  rendezvous  process  using  ICGL 
and  NG. 


(a)  Acceleration  in  the  correction  sub-phase 


Time  (s) 

(b)  Acceleration  in  the  guidance  sub-phase 

Fig.  13  Control  accelerations  in  the  correction  sub-phase  and 
guidance  sub-phase. 
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Fig.  14  UAV  and  tanker  flight  altitude  change  histories  using 
ICGL  and  NG. 


NG.  By  comparison,  it  can  be  seen  that,  during  the  correction 
sub-phase,  it  takes  52  s  for  the  UAV  to  complete  its  velocity 
direction  adjustment  using  the  ICGL  and  44.2  s  using  the 
NG.  For  the  guidance  sub-phase,  it  takes  103.4  s  to  complete 
this  process  using  the  ICGL  compared  with  200.8  s  using  the 
NG.  In  general,  it  takes  217  s  in  total  to  accomplish  the  whole 
rendezvous  process  using  the  ICGL  and  300  s  using  the  NG.  It 
is  obvious  that  the  proposed  ICGL  has  better  performance 
than  the  NG  for  the  rendezvous  problem. 

Figs.  15  and  16  show  the  change  histories  of  the  distance 
and  the  velocity  angle  between  the  UAV  and  the  tanker  during 
the  whole  rendezvous  process  using  the  ICGL,  respectively.  It 
can  be  seen  from  Figs.  15  and  16  that  the  UAV  eventually 
keeps  a  formation  flight  with  the  tanker.  This  indicates  the 
effectiveness  of  the  proposed  ICGL  method.  A b  represents 
the  distance  between  the  UAV  and  the  tanker  in  Fig.  15,  A 6 
represents  the  angle  between  the  velocity  vectors  of  the  UAV 
and  the  tanker  in  Fig.  16. 

To  test  the  robustness  of  the  proposed  ICGL  to  wind  dis¬ 
turbance,  we  add  a  step  wind  disturbance  along  the  x-axis  with 
amplitude  of  aw  =  2  m/s2  at  the  simulation  time  of  300  s  in  the 
formation  maintaining  phase  as  shown  in  Fig.  17.  The 


Fig.  15  Distance  between  UAV  and  tanker. 


Fig.  16  Angle  between  the  velocity  vectors  of  UAV  and  tanker. 


simulation  result  for  the  whole  rendezvous  process  is  shown  in 
Fig.  18.  Fig.  19  shows  the  deviations  Ax  and  Ay  of  the  UAV  to 
the  desired  distance  with  the  tanker  in  the  x-axis  and  y-axis 
directions  in  the  presence  of  wind  disturbance.  It  can  be  seen 
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Fig.  17  A  step  wind  disturbance. 
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Fig.  18  Trajectories  of  the  rendezvous  process  against  wind 
disturbance. 
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Fig.  19  Distance  between  UAV  and  tanker  on. 


that  the  steady  state  error  is  about  2  m  in  the  x-axis  direction 
and  less  than  3  m  in  the  y-axis  direction.  Fig.  20  shows  the 
change  of  the  control  acceleration  ad  of  the  UAY  under  wind 
disturbance.  The  simulation  result  indicates  that  the  proposed 
ICGL  has  good  robustness  against  wind  disturbance. 


Fig.  20  Control  acceleration  of  UAV  under  wind  disturbance. 


5.  Conclusions 

(1)  A  guidance  approach  of  the  ICGL  addressing  the  ren¬ 
dezvous  and  formation  problem  for  a  UAV  in  AAR  is 
presented  in  this  paper.  The  ICGL  divides  the  whole 
process  into  two  phases:  the  long-range  guidance  phase 
and  the  formation  phase,  in  which  the  long-range  guid¬ 
ance  phase  is  further  divided  into  two  sub-phases:  the 
correction  sub-phase  and  the  guidance  sub-phase. 

(2)  The  ICGL  solves  the  guidance  problem  of  each  phase 
and  sub-phase  in  a  two-dimensional  space  with  the  iter¬ 
ative  method  to  obtain  the  control  acceleration  for  the 
UAV. 

(3)  Simulation  results  demonstrate  that  the  proposed  ICGL 
is  effective  for  the  UAV  to  deal  with  the  rendezvous  and 
formation  problem  in  AAR  and  is  robust  against  wind 
disturbance. 


Acknowledgements 

We  would  like  to  thank  the  anonymous  reviewers  for  their 
constructive  comments  on  this  manuscript.  This  study  was  par¬ 
tially  supported  by  the  National  Natural  Science  Foundation 
of  China  (No.  61333004),  partially  by  the  Aeronautical  Science 
Foundation  of  China  (No.  20115868009),  and  partially  by  the 
open  funding  project  of  the  State  Key  Laboratory  of  Virtual 
Reality  Technology  and  Systems  at  Beihang  University  of 
China  (No.  BUAA-VR-13KF-01). 

References 

1.  Luo  DL,  Lu  TX,  Wu  SX.  A  survey  on  guidance  laws  for  flight 
vehicle.  J  Syst  Simul  2010;22(S1):  16-20  [Chinese]. 

2.  Guo  J,  Dong  XM,  Xu  YJ,  Liao  KJ.  Design  of  UAV  autonomous 
controller  for  rendezvous  in  aerial  refueling.  Contr  Decis 
20 10;25(4):  567-71  [Chinese]. 

3.  Liang  X,  Wang  HL,  Li  DW,  Lv  WT.  Three-dimensional  path 
planning  for  unmanned  aerial  vehicles  based  on  principles  of 
stream  avoiding  obstacles.  Acta  Aeronaut  Astronaut  Sin 
20 1 3  ;34(7) :  1 670-8 1  [Chinese]. 

4.  Guelman  M.  A  qualitative  study  of  proportional  navigation.  IEEE 
Trans  Aerosp  Electron  Syst  1 97 1;7(4): 337-43. 

5.  Tyan  F.  Capture  region  of  a  GIPN  guidance  law  for  missile  and 
target  with  bounded  maneuverability.  IEEE  Trans  Aerosp  Electron 
Syst  201 1  ;47(1):201— 1 3. 

6.  Lin  CL,  Li  YH.  Development  of  3-D  modified  proportional 
navigation  guidance  law  against  high-speed  targets.  IEEE  Trans 
Aerosp  Electron  Syst  20 1 3 ;49(  1 ) : 677—8 7 . 


A  guidance  law  for  UAV  autonomous  aerial  refueling  based  on  the  iterative  computation  method 


883 


7.  Zhao  SY,  Zhou  R,  Wei  C,  Ding  QX.  Design  of  time-constrained 
Guidance  Laws  via  virtual  leader  approach.  Chin  J  Aeronautics 
2010;23(1):  103-8. 

8.  Yamasaki  T,  Balakrishnan  SN,  Takano  H.  Integrated  guidance 
and  autopilot  design  for  a  chasing  UAV  via  high-order  sliding 
modes.  J  Franklin  Inst  2012;349(2):53 1—58. 

9.  Miao  SQ,  Cong  BL,  Liu  XD.  Adaptive  sliding  mode  control  of 
flexible  spacecraft  on  input  shaping.  Acta  Aeronaut  Astronaut  Sin 
20 1 3;34(8):  1 906-14  [Chinese]. 

10.  Wang  L,  Li  GC,  Wang  ZL,  Jiao  B.  Sliding  mode  control  of  an 
under  actuated  quadrotor  UAV.  J  Harbin  Eng  Univ 
20 1 2;33(  1 0) :  1 248-53 . 

11.  Jing  WX,  Li  CY,  Qi  ZG.  Application  of  the  3D  differential 
geometric  guidance  commands.  J  Astronaut  2007;28(5):  1235-40 
[Chinese]. 

12.  Peng  SC,  Pan  L,  Han  DP,  Shen  LC.  A  new  3D  guidance  law  based 
on  nonlinear  method.  Acta  Aeronaut  Astronaut  Sin 
2010;31(10):201 8— 25  [Chinese]. 

13.  Li  J,  An  HL,  Xiang  XJ,  Shen  LC.  UAVs  formation  rendezvous 
method  based  on  differential  geometry  and  Lie  group.  J  Natl  Univ 
Def  Technol  20 13;35(6):  157-64  [Chinese]. 

14.  Li  XS,  Li  X.  Robust  adaptive  backstepping  design  for  unmanned 
aerial  vehicle  formation  guidance  and  control.  J  Appl  Sci 
2012;30(5):552-8  [Chinese]. 

15.  Xia  QY,  Xu  JF.  A  design  of  triaxial  unmanned  rotor  aircraft  and 
its  adaptive  flight  control  system.  Acta  Aeronaut  Astronaut  Sin 
201 3;34(3):495— 508  [Chinese]. 

16.  Luo  YL,  Yin  CX.  Study  of  small  UAV  control  system  based  on 
fuzzy  RBF  neural  network.  J  Civil  Aviat  Flight  Univ  China 
2014;25(l):46-9  [Chinese]. 

17.  Hassan  GM,  Yahya  KM,  Haq  IU.  Leader-follower  approach 
using  full-state  linearization  via  dynamic  feedback.  Proceedings  of 
the  International  Conference  on  Emerging  Technologies',  2006  Nov 
13-14;  Peshawar,  Pakistan;  2006.  p.  297-305. 

18.  Linorman  NHM,  Liu  HHT.  Formation  UAV  flight  control  using 
virtual  structure  and  motion  synchronization.  Proceedings  of  the 
American  Control  Conference',  2008  Jun  11-13;  Seattle,  WA,  USA; 
2008.  p. 1782-7. 

19.  Karimoddini  A,  Lin  H,  Chen  BM,  Lee  TH.  Hybrid  three- 
dimensional  formation  control  for  unmanned  helicopters.  Autom¬ 
ation  2013  ;49(2) : 424-3 3 . 

20.  Han  ZZ.  An  iterative  guidance  method  for  the  large  launch 
vehicle.  J  Astronaut  1983;1:9-21  [Chinese]. 


21.  Fu  Y,  Chen  G,  Lu  BG,  Guo  JF,  Cui  NG.  A  vacuum  adaptive 
iterative  guidance  method  of  launch  vehicle  based  on  optimal 
analytical  solution.  Acta  Aeronaut  Astronaut  Sin  201 1  ;32(9): 
1696-704  [Chinese]. 

22.  Song  ZY.  From  accurate,  precise  to  perfect-manned  space 
promotes  the  development  of  guidance  method  on  launch  vehicle. 
Aerosp  Contr  20 1 3;31(1):4— 1 0  [Chinese]. 

23.  Xiao  LX,  Wang  JP,  Wang  AM,  Yang  J.  Optimal  iterative 
guidance  method  based  on  required  trajectory.  Aerosp  Shanghai 
2012;29(3):l-5  [Chinese]. 

24.  Chandler  DC,  Smith  IE.  Development  of  the  iterative  guidance 
mode  with  its  application  to  various  vehicles  and  missions.  J 
Spacecraft  Rockets  1967;4(7):898-903. 

25.  Park  S,  Deysty  J,  How  JP.  A  new  nonlinear  guidance  logic  for 
trajectory  tracking.  Proceedings  of  the  AIAA  Guidance,  Navigation, 
and  Control  Conference  and  Exhibit',  2004  Aug  16-19;  Providence, 
Rhode  Island,  USA;  2004.  p.  16-9. 

26.  Luo  DL,  Li  YF,  Wu  WH,  Shen  CL.  Terminal  guidance  law  based 
on  receding  horizon  control  for  homing  missile.  J  Nanjing  Univ 
Aeronaut  Astronaut  2005;37(l):52-6  [Chinese]. 

Luo  Delin  received  Ph.D.  degree  in  navigation,  guidance  and  control 
from  Nanjing  Universtiy  of  Aeronautics  and  Astronautics  in  2006.  He 
is  presently  an  associate  professor  with  the  Department  of  Automa¬ 
tion,  School  of  Information  Science  and  Technology,  Xiamen  Uni¬ 
versity.  His  research  interests  include  flight  vehicle  guidance  and 
control,  UAV  cooperative  control  and  decision  making. 

Xie  Rongzeng  received  B.S.  degree  in  automation  from  Wuhan  Uni¬ 
versity  of  Technology  in  2010.  Currently,  he  is  a  M.S  candidate  at  the 
Department  of  Automation,  Xiamen  University.  His  research  interests 
include  UAV  guidance  and  cooperative  control. 

Duan  Haibin  received  Ph.D.  degree  in  control  theory  and  control 
engineering  from  Nanjing  Universtiy  of  Aeronautics  and  Astronautics 
in  2005.  He  is  currently  a  professor  with  the  School  of  Automation 
Science  and  Electrical  Engineering,  Beihang  University.  His  research 
interest  include  flight  control,  autonomous  control  of  multiple  aerial 
vehicles,  computational  intelligence. 


